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The superconducting state emerges due to the formation and condensation 
of carrier Cooper pairs, although the exact microscopic mechanism responsi- 
ble for the pairing in different materials varies and, in many cases, remains 
elusive. An important step towards a microscopic understanding is detailed 
knowledge of the superconducting order parameter. This is especially impor- 
tant in materials such as Sr 2 Ru0 4 (1, 2) and UPt 3 (3) where the carriers in 
the Cooper pairs are believed to form a triplet, introducing an additional de- 
gree of freedom and complexity to the problem. As a result the exact nature 
of the order parameter in these compounds, which have become paradigms 
for unconventional superconductivity, is still unresolved. We have studied the 
vortex lattice in Sr 2 Ru0 4 to measure the intrinsic anisotropy of the supercon- 
ducting state between the c axis and the Ru-O basal plane (~ 60), and found 
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this to greatly exceed the anisotropy of the upper critical field (~ 20). Our 
result imposes significant constraints on possible models of triplet pairing in 
Sr 2 Ru0 4 . 

Introduction Multiple experimental and theoretical studies provide compelling support for 
triplet pairing of carriers (electrons and/or holes) and an odd-parity, p-wave order parameter 
symmetry in superconducting Sr 2 Ru0 4 . However, seemingly contradictory experimental results 
have left important unresolved questions concerning the detailed structure and coupling of the 
orbital and spin parts of the order parameter. A prominent example of this predicament is 
conflicting evidence as to whether the p-wave order parameter is chiral (4, 5). 

In this Report we focus on the anisotropy of the superconducting state of Sr 2 Ru0 4 . The 
Fermi surface in this material consists of three largely 2-dimensional sheets with Fermi velocity 
anisotropics ranging from 57 to 174 (1, 6), and one would thus expect an upper critical field 
anisotropy within this range (7, 8). Experiments, however, find a much smaller upper critical 
field (H c2 ) anisotropy T Hc2 = H^f/H^ — 20 and a near constant H c2 when the applied field is 
within ±2° of the basal plane (9). Within the same angular range the superconducting transition 
at H c2 becomes first-order, leading to suggestions of a subtle coupling between the magnetic 
field and the triplet order parameter (10), or Pauli limiting which is inconsistent with a Cooper 
pair zero spin-projection along the c axis (11). 

We used small-angle neutron scattering (SANS) studies of the vortex lattice (VL) to deter- 
mine the intrinsic anisotropy (r ac ) of the superconducting state in Sr 2 Ru0 4 yielding T ac ps 60. 
This is three times greater than V Hr2 , and poses significant constraints on the possible order 
parameter symmetry in this compound. 



Materials and Methods The experiment was performed using a single crystal of Sr 2 Ru0 4 
grown by the floating zone method and carefully annealed, yielding a critical temperature T c = 
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1.45 K and no indication of a 3 K-phase (1). Measurements were performed at T = 40 — 60 mK 
using a dilution refrigerator (DR) inserted into a horizontal-field cryomagnet. Magnetic fields 
of HqH = 0.5 T and 0.7 T were applied close to the sample a axis. A motorized fi-stage could 
rotate the DR within the magnet, allowing in-situ sample alignment and measurements as the 
crystalline basal plane was rotated with respect to H. A schematic of the experimental configu- 
ration is shown in Fig. 1(A). The VL was prepared by changing H and Vl at the base temperature, 
followed by a damped small-amplitude field modulation. This method produces a well-ordered 
VL and eliminates the need for a field-cooling procedure before each measurement. The SANS 
experiment was carried out on the Dll and D22 instruments at Institut Laue-Langevin (ILL), 
using a neutron wavelength A n = 1.7 nm and a wavelength spread AA„/A n = 10%. Some mea- 
surements were performed using polarized incident neutrons and a 3 He analysis cell to allow 
discrimination between spin-flip and non spin-flip scattering. 

Results The VL scattered intensity in a SANS experiment is determined by the amplitude 
of the field modulation and is proportional to |h| 2 , where h(q) is the Fourier transform of the 
magnetic field B(r) (12). Using state-of-the-art SANS instruments at a high-flux neutron source 
such as ILL, it is possible to measure the diffraction from a well-ordered VL with a longitudi- 
nal Fourier coefficient \h z \ as low as 0.1 — 1 mT, depending on the amount of background 
scattering (13). Here \h z \ oc Aj 2 , where Aj_ is the average penetration depth in the screening 
current plane perpendicular to the applied field. Previous SANS studies with H || c found a 
VL form factor for Sr 2 Ru0 4 no greater than a few mT (14). This indicates that measurements 
with H 1 c should not be possible as \h^ c \ 2 /\hf\ 2 oc (A a6 /A c ) 2 = T" 2 , and with T ac > 20 
we estimate \h^ c \ < 3 fiT, at least two orders of magnitude below what is required for a VL 
SANS experiment. However, in highly anisotropic superconductors such as Sr 2 Ru0 4 , there is a 
strong preference for the vortex screening currents to run within the basal afr-plane. As a result, 
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a small angle f2 between the applied field and the basal plane will lead to a significant trans- 
verse Fourier coefficient (h x ). Estimates based on an extended London model which includes 
an effective mass anisotropy yields \h x /h z \ 2 oc T 2 ac (15), and thus predict h^ c to be comparable 
in magnitude to h\ c . As a result, scattering due to the transverse field modulation should be 
observable. This is confirmed by the VL diffraction pattern shown in Fig. 1(B) which show 
Bragg peaks aligned with the crystalline b direction (y axis), is in contrast to the exotic VL 
configurations predicted by Agterberg (16) and Kita (17). 

Scattering from the transverse field modulation leads to a flip of the neutron spin (a _L h x ) 
and a Zeeman- splitting of the VL rocking curves shown in Fig. 2(A). Two maxima are observed 
for both the top (positive Q y in Fig. 1(B)) and bottom (negative Q y ) VL reflection, as the angle 
(</?) between the scattering vector Q and the direction of the incident neutron beam is varied to 
satisfy the Bragg condition (18). As expected, no scattering from the otherwise more commonly 
observed longitudinal VL field modulation (h z ) could be measured in Sr 2 Ru04. 

To verify that the observed diffraction is due to spin-flip scattering, measurements with a 
polarized neutron beam were performed. Here only one maximum is observed for each Bragg 
reflection, selected according to the direction of the neutron spin (18). Furthermore, the scat- 
tered intensity normalized to the incident neutron flux is doubled relative to the unpolarized 
beam as expected. Moreover, using polarization analysis it is possible to measure only the 
spin-flip scattering as shown in Fig. 1(B). 

Dividing the integrated intensity by the incident neutron flux yields the integrated VL re- 
flectivity 

where 7 = 1.913 is the neutron magnetic moment in nuclear magnetons, t is the sample thick- 
ness and $0 — h/2e — 2069 Tnm 2 is the flux quantum (13). As shown in Fig. 2 each peak is 
fitted to the sum of three Gaussians due to the asymmetry of the rocking curves (19). Moreover, 



the integrated intensity for the two maxima (top, bottom) for a given reflection are added, as 
each corresponds to half the incident flux (one direction of the neutron spin). The form factor 
obtained in this fashion is shown in Fig. 3, for all measured fields and JTs. This illustrates how 
the VL SANS measurements are possible within a narrow angular range, with H close to, but 
not perfectly aligned with, the basal plane. The width of the measurement "window" decreases 
with increasing field due to the rapidly decreasing H c2 (VL) (9). In addition, the overall form fac- 
tor decreases with increasing field. While the anisotropic London model provides a qualitative 
description of the enhanced field modulation (15), it is not able to provide a good quantitative 
fit to the data. As shown in Fig. 3, an extended London model that include a so-called core 
correction by multiplying the calculated \h x \ by exp(— cQ 2 (Vt) does still not yield a good 
fit to the data. Here c ~ 1, Q(Cl) is the magnitude of the VL scattering vector (see below) and 
£,ab = (^o/ZkH^) 1 / 2 is the in-plane coherence length (13). A quantitatively accurate model for 
the VL form factor is highly desirable as it would allow a determination of both A and £. 

We now turn to the main result of this Report, which is the measurement of the VL anisotropy. 
In an anisotropic superconductor the VL Bragg peaks are expected to lie on an ellipse with a 
major-to-minor ratio given by (7) 



as shown in Fig. 4(A). This ^-dependence was derived for anisotropic (but still 3 -dimensional) 
superconductors, and was verified in early VL SANS measurements on 2ff-NbSe 2 with r ac = 

3.2 (20). Although Sr 2 Ru04 is a layered material, the coherence length along the c axis £ c = 

3.3 nm is still several times greater than the Ru-0 interlayer spacing (0.64 nm) (1), and we 
expect Eq. © to be applicable (21). 

Due to the large anisotropy in Sr 2 RuOzi, VL Bragg peaks which are not on the vertical axis 
have scattering vectors essentially parallel to h x making them unmeasurable as only compo- 



r, 



ac 




(2) 



5 



nents of the magnetization perpendicular to Q will give rise to scattering (22). Instead we 
determine the VL anisotropy based on flux quantization. Assuming that each vortex carries one 
flux quantum $ , the area of the reciprocal space ellipse in Fig. 4(A) is determined uniquely 
by the applied magnetic field. This yields Y VL = (Qo/Q) 2 , where Q is the magnitude of the 
measured VL scattering vector and Q = 2ir(2fi H / y/S&o) 1 / 2 corresponds to an undistorted 
hexagonal VL (r ac = 1). The magnitude of Q can be determined either from the position of the 
VL Bragg peaks on the detector as shown in Fig. 1(B), or from the peak positions <p\, . . . , (p± in 
Fig. 2(A) (18). The two methods yield nearly identical results and using the average Q we ob- 
tain T VL (n) shown in Fig. 4(B). Within the scatter in the data the results for both fields collapse 
onto a single curve, increasing upon approaching the a axis and reaching a value slightly above 
50 before the intensity vanishes. If one assumes a quantization of $o/2, as recently reported 
for mesoscopic rings of Sr 2 Ru0 4 (25), the deduced values for T VL would double. However, 
we consider this an unrealistic scenario in the present case, with a macroscopic, uniform sam- 
ple. Fitting the data to Eq. © yields r ac = 58.5 ± 2.3. Only for angles within ±1.3° does 
the measured anisotropy deviate from that expected for an infinite ac anisotropy. Also shown 
for comparison is Tvl expected from the low temperature r# c2 = 20 (9). We note that r ac 
coincides with the anisotropy of the j3 Fermi surface sheet (57) (1, 6). 

Discussion The large difference between r# c2 and the intrinsic anisotropy of the supercon- 
ducting state deep within the mixed phase measured by r ac , indicates a strong suppression of 
the upper critical field in Sr 2 Ru04 for H _L c. One possible explanation for this difference is 
Pauli limiting due to the Zeeman splitting of spin-up and spin-down carrier states by the applied 
magnetic field and the resulting reduction of the superconducting condensation energy (24). In 
spin-triplet superconductors the order parameter is most conveniently described in terms of the 
d- vector, directed along the zero spin projection axis where the configuration of the Cooper 
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pairs is given by l/v^dtl) + lit)) (1> 2, 5). Consequently, Pauli limiting in the triplet case can 
only occur when H || d. This is inconsistent with the chiral superconducting state proposed for 
Sr 2 Ru0 4 which requires d || c (2, 5), at least for the current experimental configuration with 
H JL c. 

Instead our results indicate that the superconducting state in this material should be un- 
derstood based on either a subtle coupling between the magnetic field and the triplet order 
parameter (10); recent multiband p-wave models (25); a field-dependent mixing of singlet and 
triplet states (26) or singlet superconductivity (11). It should be noted however, that the latter 
(singlet pairing) does not provide an satisfactory explanation for the extreme sensitivity of T c to 
impurities or to the chiral properties of Sr 2 Ru0 4 (1, 2). 

Conclusion In conclusion, we have used SANS to measure the anisotropy of the supercon- 
ducting state in Sr 2 Ru0 4 which greatly exceeds that of the upper critical field. This result 
imposes significant constraints on the possible pairing of carriers in this material. Any model 
aimed at describing the superconducting phase must provide a satisfactory explanation for this 
observation as well as the temperature dependence of T Hc2 , which increases with T and reaches 
T ac at T c (27), and the first-order nature of H c2 close to the basal plane (10). 
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Figure 1 : (A) Experimental geometry. The coordinate system is defined with z along H and y 
in the Ru-0 basal plane (along b). The applied magnetic field H is rotated away from the Ru-0 
(spanned by a and b) by an angle f2. Neutron spins (a) are are parallel or antiparallel to the 
magnetic field. The incident neutron beam is in the yz-plane, at an angle <p relative to the field 
direction. The observed VL scattering vector is denoted Q and the longitudinal and transverse 
component of the field modulation by h z and h x respectively. (B) Diffraction pattern showing 
spin flip scattering from the VL due to the transverse field modulation (h x ). The two Bragg 
peaks corresponds to ±Q in the panel above. No background subtraction was performed but a 
small remnant of the undiffracted beam (due to the finite flipping ratio ~ 8) close to Q = is 
masked off. 
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Figure 2: Vortex lattice rocking curves showing the scattered intensity as a function of the angle 
cp, for an unpolarized neutron beam (A) and for neutrons polarized with their magnetic moment 
parallel (B) and antiparallel (C) to the applied field. Except where shown, error bars are equal 
to or smaller than the symbols. With unpolarized neutrons two maxima are observed for both 
the bottom (^1/3) and top ((^2/4) VL Bragg reflections. In all cases the intensity was normalized 
to the incident neutron flux. The curves are fits to the data as described in the text. 
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Figure 3: Vortex lattice form factor at 40 mK as a function of applied field and angle with the 
a axis. The statistical error is roughly the size of the symbols. The solid lines are guides to the 
eye. The dashed line show an extended London model fit to the 0.5 T data as discussed in the 
text, with X ab = 650 nm, £ afe = 66 nm, c = 1/2 and r ac = 58.5. 
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Figure 4: Vortex lattice anisotropy. (A) Schematic of VL Bragg reflections lying on an ellipse 
with major-to-minor axis ratio, Tvl = 6. Only the filled (red) peaks are observed. The recipro- 
cal space area of the ellipse is nQl = 87i 3 fi H/ a/3$o- (B) Measured VL anisotropy at 40 mK 
as a function of applied field and angle with the a axis Except where shown explicitly the 
statistical error is the size of the symbols. The lines show the VL anisotropy calculated using 
Eq. © and T ac = 20 (dotted), 58.5 (dashed) and oo (full). 
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The two different directions of the neutron spin with respect to the applied field correspond 
to different nuclear Zeeman energies and lead to opposite shifts of the neutron momentum 
vector 



k Ui) = k o V 1 ± Ae / £ o, (!) 

where the subscript in parentheses henceforth correspond to the lower (minus) sign in the ± Ae 
term. Here the nominal neutron wavevector k = 2n/\ n , Eq = h 2 kQ/2m n , Ae = j/^nB and 
7 = 1.913 is the neutron magnetic moment in nuclear magnetons /xjv = eh/2m n = 31.5 neV/T. 
With A„ = 1.7 nm and /i Q H = 0.5 T one finds jfe| - jfej = 2 x 10" 4 fcg. 

Due to the short Q in the range 0.003 — 0.01 k , the small shift in the neutron wavevector 
nonetheless leads to a significant difference in the angle (<£>i(2)) required to satisfy the Bragg 
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condition as shown in Fig. S2. In this case Bragg's law is replaced by 

fcf - k\ ± Q 2 = 2k m Q sin ip 1{2) . (2) 

In contrast the scattering angles 9i/8 2 = k^/k^ are essentially identical for the two scattering 
configurations. 

The magnitude of Q can be determined from the peak positions (pi, . . . , </? 4 (see Fig. 2(A) in 
the main text) by 

Qi(2) = =F^ta) sin <Pi(2) =F ~ k t(i) cos 2 ip 1{2) (3) 

Qs(4) = Th(V Sil1 ^3(4) ± V / ^tO) ~ k W) COs2 ^3(4), (4) 

using the same notation as earlier. 

A more detailed discussion of the spin-flip scattering can be found in Ref. (1), where a 
similar but less extreme effect was observed in YBCO. 
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Figure SI: Schematics showing the scattering geometries corresponding to the reflection at 
cp — (fx (A) and (p 2 (B). Not to scale. 
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